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ABSTRACT 

We argue that it is currently impossible to simulate X-ray clusters using correct 
equations, because even the MHD description is not applicable. But since fluid 
simulations actually reproduce observations quite well, one may try to improve the 
fluid codes by including molecular transport of heat and momentum. We calculate the 
effective molecular viscosity for the simplest model of magnetic field and obtain 1/5 of 
the Braginskii value, similar to 1/3 of Spitzer for the heat conduction. This is large 
enough to noticeably damp the X-ray cluster turbulence. 



1. Introduction 

To turn X-ray clusters into a precise probe of cosmology, one needs to simulate the clusters 
numerically. This is being done (eg. Nagai, Kravtsov, Vikhlinin 2006, Sijacki, Springel 2006 and 
references therein) but there are some fundamental problems with these or any potentially doable 
simulations (§2). 

Given that the first principle numerical simulation is not possible, one can still try to use 
some effective fiuid description. And in fact remarkable agreement with observations was already 
achieved (Nagai, Kravtsov, Vikhlinin 2006). Further improvement, which is needed for cosmology, 
may come when the missing physics is included. One obvious thing is to include molecular heat 
conduction and viscosity (eg. Gruzinov 2002, Sijacki, Springel 2006). 

It has been shown that in the simplest model with random magnetic fields, the effective 
molecular heat conduction equals 1 /3 of the unmagnetized (Spitzer) value. Here we show (§3) that 
the same magnetic field model gives effective molecular viscosity equal to 1 /5 of the unmagnetized 
(Braginskii) value. Of course, practical numerical simulations should be done using some fudge 
factors correcting the Braginskii-Spitzer transport coefficients (Sijacki, Springel 2006). We are 
merely pointing out that 1/3 and 1/5 are, in principle, realizable values of the fudge factors. 



2. Fluid description of X-ray clusters 



Fluid description is not applicable to the cluster plasma. Here by the fluid description we 
mean ideal or non-ideal, pure hydro or MHD, or even Braginskii (1965) non-ideal multiple-fluid 
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MHD. Collisional Vlasov (Maxwell electrodynamics plus Boltzmann equation with collisions) is, 
of course, valid, but it requires particle simulations which seem to be well out of reach. 

The problem with the fluid description is the magnetic field. Even if it turns out that the 

magnetic field is very weak (and this is a big if given the desired ^1% accuracy), one must still 
calculate it, because the magnetic field determines the plasma transport coefficients. And it is 
known that molecular heat conduction and viscosity, if close to the Braginskii-Spitzer values, have 
a significant effect on the cluster (eg. Gruzinov 2002, Sijacki, Springel 2006). 

If magnetic fields must be calculated, why cannot one calculate those using the Braginskii 
equations? This is best explained by the frozen-in law. 

According to Maxwell equations (Faraday law), S^B = V x E. For the cluster plasma, 
the electric field component parallel to the magnetic field must nearly vanish (plasma electric 
conductivity along the magnetic field is very high). But if E is perpendicular to B, one can 
introduce an auxiliary vector field v, such that E = v x B, and then write the Faraday law as 



This equation means that the magnetic field B is frozen into the auxiliary field v. Due to the 
frozen-in law, the magnetic field topology (in particular the magnetic helicity / d^rA. ■ B, where A 
is the vector potential, B = V x A) is conserved. 

If one could calculate the auxiliary field v (only the components of v perpendicular to B are 
needed and defined), the frozen-in equation (1) would give the magnetic field. On large scales, 
larger then the plasma mean free path 



the auxiliary field v coincides with the plasma velocity. The problem is that the magnetic field 
still exists and evolves iso-topologically on length scales of order A and smaller. For example, on 
the length scale equal to the plasma mean free path, the magnetic field damping time tmd (due to 
finite electric conductivity) can be written as 



where Tg is the electron collision time, and r^, is the classical electron radius. This number is too 
big to be taken literally. Kinetic plasma instabilities will probably develop and change the field - 
but the problem is that we don't know exactly how the field will evolve on these length scales. 

We can also state the problem with the fluid description in a slightly different language. 
As a warm-up, consider fluid description of an unmagnetized gas. The gas is characterized by 
the mean free path A and the sound speed Cg. Suppose a flow of velocity ^ V < Cg is imposed 
at length scales ~ L » A. According to the fluid description, Kolmogorov cascade to smaller 



5tB = V X (v x B). 



(1) 




(2) 




(3) 
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scales develops, with characteristic velocity vi ~ y(Z/L)^/^ on all scales I from the inertial 
range < I < L. Here the dissipation scale is given by the equality of the eddy frequency 
and the damping rate: z^/Z^ ~ vi^/l^i, where u ~ c^X is the molecular viscosity. It follows that 
Id/X ~ {cs/V)^^^{L/\)^^^ ^ 1. This means that the fluid description is consistent - the flow is 
fully damped on scales much larger than the mean free path A. 

Now consider magnetic fields in the cluster plasma, assuming Braginskii MHD description. We 
don't know how to calculate magnetic field cascade to small scales. We will simply underestimate 
the characteristic "eddy" frequency replacing it with the large scale frequency V/L. Then the 
magnetic dissipation scale Imd is given by i^m/lmd ~ ^/-^j where u^n ~ c^/a is the magnetic 
diffusivity, and a ~ c^reUTe is the electric conductivity. This gives 

Again this is too small - kinetic effects will turn on and cut off the cascade. But we don't know 
what actually happens. 

The problem of the cluster magnetic fields and molecular transport is related to the well-know 
problem of high-energy astrophysics - effective magnetohydrodynamical description of collisionless 
plasma. It might turn out that MHD (with anisotropic molecular transport) can mimic the 
outcome of kinetic instabilities, but this has not been demonstrated. 



3. Effective viscosity in a random time-varying magnetic field 

We start by repeating our (trivial) argument for the Spitzcr over 3 heat conduction. When 
the same procedure is applied to viscosity, one gets Braginskii over 5. 

Heat conduction: Heat conduction along the magnetic field is equal to KsVyT, were VyT 
is the temperature gradient along the field, ks is the Spitzer thermal conductivity (Braginskii 
1965). If the magnetic field evolves randomly and fast ^, one should average over the magnetic 
field direction. Averaging the heat flux over the random magnetic field direction b gives Spitzer 
over 3: 

<q>= - <b/«5b- Vr>= -^VT. (5) 

Viscosity: Now we should average the momentum flux tensor 

3 1 

T^ij = -VB-^ihbj - -Sij)bnbmWnm, (6) 



Say, the cluster plasma is being crisscrossed by about a hundred galaxies. The galaxies which have retained their 
plasma entrain the magnetic field, forcing reconnection on the short - galaxy crossing - time scale. 
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where the rate of strain tensor is 

Wn^ = dnVm + 5„VW " • Y5nm. (7) 

and V is the plasma (ion) velocity. Averaging eq. (6) over the random unit vector b gives Braginskii 
over 5: 

< T^ij >= -^Wij (8) 

We now re-derive the suppression factors 1/3 and 1/5 using the heat production rate. For 
heat conduction, Q oc (VT)^ without magnetic field, and Q oc {h ■ VT)^ with magnetic field. 
Averaging over b gives the factor 3 suppression. 

For viscosity, Q = (775/2) (W)^, without magnetic field. With non-zero magnetic field, the 
heat production rate must be proportional to (b • W ■ b)^, because all other components of the 
rate of strain can be represented by the velocity field with a gradient perpendicular to b, giving a 
suppressed transport and damping. Then Q = (?7B/2)(3/2)(b • W • b)^. Here the normalization 
coefficient 3/2 can be obtained by taking the rate of strain tensor W = diag(— 1/2, —1/2, 1) for 
b = (0, 0, 1), and noting that for this geometry there should be no suppression of heating (without 
collisions, parallel and perpendicular temperatures change in the same way with and without 
magnetic field). Averaging the magnetized heating rate over b gives the factor 5 suppression. 

I thank Alexey Vikhlinin for discussions. This work was supported by the David and Lucile 
Packard foundation. 
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